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We present a proposal for the estimation of Biittiker-Landauer traversal time based on the visibil- 
ity of transmission current. We analyze the tunneling phenomena with a time-dependent potential 
and obtain the time-dependent transmission current. We found that the visibility was directly con- 
nected to the traversal time. Furthermore, this result is valid not only for rectangular potential 
barrier but also for general form of potential to which the WKB approximation is applicable . We 
compared these results with the numerical values obtained from the simulation of Nelson's quan- 
tum mechanics. Both of them fit together and it shows our method is very effective to measure 
experimentally the traversal time. 

PACS number(s): 03.65.Bz, 73.40. Gk 



I. INTRODUCTION 

Soon after the advent of quantum mechanics, MacCoU suggested that there is a time associated with the passage of 
a particle under a tunneling barrier, i.e. a tunneling time [1]. Now the time has been measured in several experiments 
and its qualitative results have been obtained. However, it is not clear whether a unique time exists or not, since we 
have no univocal definition of tunneling time and no definite experimental data. See [2] , [3] and references therein for 
reviews of the problem. 

In this paper, we present a proposal for the estimation of Biittiker-Landauer traversal time based on the visibility 
of transmission current experimentally. Biittiker and Landauer [4], [5] invoked an oscillatory barrier to estimate a 
tunneling time. The original static barrier was augmented by a small oscillation in the barrier height. The amplitude 
of the oscillation is kept small; the disturbance of the original kinetics can be made small as desired. At very low 
modulation frequencies the incident particle sees a particular part of the modulation cycle. The particle sees an 
effectively static barrier, but later parts of the incident wave see a slightly different barrier height. As one turns up 
the modulation frequency, one eventually reaches a range where an incident particle no longer sees a particular portion 
of the modulation cycle, but is affected by a substantial part of the modulation cycle, or several cycles. They claimed 
that the frequency at which this transition occurs, i.e., the frequency where one begins to deviate substantially from 
the adiabatic approximation, is an indication of the length of time that a particle interacts with the barrier. They 
made carefully several comments as follows: It is, of course, an approximate indication of a time scale. It is not the 
eigenvalue of a Hamiltonian, indicative of a precisely measurable value. Moreover, this traversal time value may really 
be characteristic of a statistical distribution. 

They showed that for an opaque rectangular barrier, the modulated barrier approach yields 

T = dm/hn , (1) 

where d is the barrier length and Hk the magnitude of the imaginary momentum under the barrier. For a potential 
that allows the WKB approximation, it yields 

where B means the barrier region. 
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This gives a plausible estimation of traversal time based on a theoretical background. However, if one wants to 
measure the value of traversal time by an experiment, one has to draw it from the asymptotic behavior of transmission 
rate as a function of lu. Generally its dependence on to does not change so rapidly, that one cannot easily estimate 
the value from experimental data. There is another type of experiment; one projects a stationary incident particle 
beam on the target with oscillating barrier and measure the time dependence of transmission current which may also 
oscillate with the same frequency. Here we show the visibility of oscillating current gives us a good information about 
traversal time. 



II. TIME-DEPENDENT BARRIER 

Following [4] , [5] and [6] , we start by considering a Hamiltonian, 

where Va{x) is static and Vi{x) is the amplitude of a small modulation. Incident partickjs with energy E interact- 
ing with the perturbation Vi cos cut, will emit or absorb modulation quanta hu). The Schrodinger equation of this 
Hamiltonian has the solution in the barrier region 

/ rplf\ "=oo / V \ 

nx,t;E')=^E'{x)e^p(-i—] J2 ^" ( ^ ) (4) 

^ ^ n— — oo ^ ^ 

where 4>e'{x) is an eigenfunction of the time-independent Hamiltonian Hq = —{h^/2m)(f/dx'^ + Vq, HQ(j)E' = E'(j)Ei 
and Jn is a Bessel function. The time modulation of the potential gives rise to sidebands describing particles which 
have absorbed (n > 0) or emitted (n < 0) modulation quanta. Therefore we have to take into account the many 
sidebands of which the Bessel functions are appreciable. 

To the left of the barrier, we allow an incident wave at energy E and reflected waves at energies E' = En = E+nhw, 

*i(x,t) =e*'=^e-^**+ AnQ-'^-'-'e-'^', (5) 



where fc„ = y , Eq = E and fco = k. See Fig.l. We consider only the positive energy solutions. In the barrier 

region, in addition to the solution (4) with E' = E, there exist other evanescent (and oscillating, in a certain case) 
modes corresponding to the reflected wave with energy E' = En- Here we also consider only positive energy solutions. 
Taking account of these points, we have a solution in the barrier region. 



"ctt 



where k„ = xl ^"^(^^ . For the transmitted wave, we have 



*™(a;,i)= Dne'''"^e-'^K (7) 

E„>0 

For small Vi, Jn is proportional to (Vi/2?!cj)" and thus, only the small numbers of terms in the summation of (6) 
contribute effectively. Correspondingly the numbers of terms in the summations of (5) and (7) are suppressed. To 
find the solution for the Schrodinger equation, we match a superposition of incident and reflected waves (5), and also 
transmitted waves (7), at each energy En, to solutions within the barber (6). As a result of somewhat tedious but 
straight calculation(see the Appendix A), we have the transmission and reflection coefficients in the leading order, 

^ _ Jn{Vi/hLo) 2Doe'(>'->'-W2 
JoiVi/hoj) det{kn,Kn) 

X {{k^ — knko) sinh Knd — {k^ — knko){Kn/Ko) sinhnod 

+ iKn{kn + fco)(coshK;orf— coshK„d)} , (8) 
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and 

" Jo{Vi/hu)) det(fc„,K„) 

X {{k^ — knko) sinh Knd cosh Kod 
— {k^ + knko){Kn/ ko) cosh Knd sinh Kod 
+iKn{ko — kn){'i- — cosh Knd cosh Kod) 

—i{{koKn/Ko) — knKo) sinh Knd sinh Kod} , (9) 

^Kn ^" ikji^C {Kn ^'^n) 

= 2(k^ — kn) sinh Knd — AiknKn cosh «;„rf. (10) 

From these results we can obtain the transmission probability defined by the ratio of transmitted current jm and the 

incident current jinc = hk/m. It depends on the time as well as the position of measurement due to the interference 
among different energies waves. However, if we take a time average of the ratio, its dependence will disappear, 

"ett 7 

n=0 

We show an example of numerical result of the time-averaged transmission probability in Fig. 2. 

Now we will discuss the traversal time. As following to Biittiker and Landauer, we assume that Uw <C E, so that 
the wave numbers of the sidebands are 



where det(fc„,K;„) is defined by 

det(A!yi, Kn) 



fc,„ = V ^-(^^^M ^fc±n^, (12) 



and assume Uw <^Vo — E, so that 



2m{Vo-ETTy^) rriLO 

In the case of opaque barrier, taking account of the asymptotic forms of transmitted wave amplitudes, 

£•±1 = ±^£'o(e±'^^ - l)e=F^^ , (14) 

Biittiker and Landauer included first order corrections to the static barrier and obtained the intensity for the trans- 
mitted sidebands, for the case of small Vi , 

where r = md/fiK. From this expression they found that there exists the crossover from the low frequency behavior 



k±i fViT 
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where the two intensities of the sidebands are equal, to the high frequency behavior 



+^ ko \2huj J 



To, (17) 



%imVo. (18) 



ko \ 2hu) 

where the two intensities differ strongly. This transition to imbalance is best described by 



fc-iT+i + fc+iT_i 



tanhwr. (19) 



Thus they claimed the crossover from the low frequency behavior to the high frequency behavior yields the traversal 
time. 
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III. VISIBILITY AND TRAVERSAL TIME 



Their claim is a very interesting idea to estimate a certain kind of tunneling time, but it is rather difficult to 

determine its value from experiments. Now let us consider the time dependence of the transmitted currents. If one 
observes the currents at a fixed point x = L, one may see the interference effect between the different frequency waves 
in the first order approximation, 



i{k-iL-E-it) 



1 



IDor + ^l-Dol {{ko + fc+i)|-D+i| + (fco + k-i)\D_i\) cos(a;t - 

/Co 



(20) 



where (j>{L) is a phase which is independent on time t, 



<j){L) = (j)+i{L) = arg 



^0 



= -(I)--l{L) = - arg 



+ (fc+i - kQ)L 



Do 



+ (fc_i - ko)L 



(21) 



Here the asymptotic forms (14) were used. Now we show the numerical result of the time dependence of transmitted 
current at a fixed point in Fig.3. If a detector has a good time resolution, one may measure this visibility of the 
transmitted wave 



-^vis — 



max J- mm 



max t ' mm 



= ^((^o + fc+i) 



D 



+1 



+ {ko + k-i) 



D. 



Dn 



In the case of a small perturbation Vi and an opaque static potential, equation (22) is approximated by 

/v 



2Fi . , 
— — smh CUT, 



from which the traversal time is expressed by the visibility as follows; 



1 



T = — sinh 



flLO 



If one can choose an experimental setup satisfying the condition wt ^1, this expression becomes to 

2^1 ™- 



(22) 



(23) 



(24) 



(25) 



For the case of general potential shown in Fig. 4, which allows the WKB approximation, we have a transmitting 
wave after the potential wall, 



~iEt/n 



1 + 



Ji (Vi/huj) kojx) 

Jo{Vi/rioj)V ki{x) 



exp < I 



?iko{x' 



-dx' \ (1 - Si)e 



+ 



J_l (Vl/huj) ko{x) 



JoiVi/hu;) V k-iix) 



cxp< -I 



hko{x') 



-dx' \ (1 - S_i)e'' 



(26) 



where 
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4 + S±i 



Sn = exp ^- K;„(x)(ia;^ . (28) 

The detailed calculation is given in Appendix B. For an opaque potential, the damping factors 5„ are so small, that 
the transmitted current becomes 

T^S',i^l + 2jM^/^(S_i - El) cos M - 0(x))| , (29) 



where 



Therefore the visibility is given by 



_ JliVl/hu^) 



-i2sinh — / -^dx], 31 



and Eq. (25) is replaced by the following expression, 

2Vi 



m f^^ 1 h 
-TwKB = ^ / — r^dx ~ -^^^Ivis- (32) 



IV. COMPARISON OF NUMERICAL RESULTS WITH THE SIMULATION 
BASED ON THE NELSON'S QUANTUM MECHANICS 

Here we evaluate the tunneling time by the use of Nelson's approach of quantum mechanics [7] and compare them 
with numerical results of traversal time obtained from the visibility. Nelson's quantum mechanics, using the real-time 

stochastic process, enables us to describe individual experimental runs of a quantum system in terminology of the 
"analog" of classical mechanics, i.e., the ensemble of sample paths. These sample paths are generated by the stochastic 
process, 

dx{t) = iu{x{t), t) + v{x{t), t))dt + dw{t), (33) 

where x{t) is a stochastic variable corresponding to the coordinate of the particle, and u{x{t),t) and v{x{t),t) are 
the osmotic velocity and the current velocity, respectively. The dw{t) is the Gaussian white noise with the statistical 
properties of 

{dw{t)) = 0, and {dw{t)dw{t)) = —dt. (34) 

m 

In principle the osmotic and the current velocities arc given by solving coupled two equations, i.e., the kinetic equation 
and the "Newton- Nelson equation". The whole ensemble of sample paths gives us the same results as quantum 
mechanics in the ordinary approach. Once the equivalence of Nelson's framework and ordinary quantum mechanics 
is proved, it is convenient to use the relation 

h d h d 

u = Re — —lnih(x,t), and v = lm — —IntbixA), (35) 
m ox m ox 

where tp is the solution of Schrodinger equation. Since individual sample path has its own history, we obtain infor- 
mation on the time parameter, e.g., the traversal time [8], [9]. 

Now using the Nelson's quantum mechanics, we estimate the traversal time crossing over a time-dependent potential 
barrier shown in Fig.l. Suppose a simulation of tunneling phenomena based on (33), starting from t = — oo and ending 
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at t = oo. As wc treat a wave paeket satisfying the time-dependent Schrodingcr equation, the wave packet is located 
in region I initially and turns finally into two spatially separated wave packets which are in regions I and III. Fig. 5 
shows a typical transmission sample path calculated by Eq.(33) with "backward time evolution method" [8], [9]. The 
traversal time using this approach, TNeison, is defined as the averaged time interval in which the random variable x{t) 
stays in the barrier region II. Thus TNeison defined in this way has a character of statistical distribution as pointed in 
[4], [5], since it is the value averaged over the ensemble of sample paths having the transmitting wave packets. 

We call the traversal time obtained by the visibility of transmission current, Tvis- Let us compare Tvis with TNeison 
and TwKB in a rectangular potential barrier numerically. Here we take the unit with m = Ti = 1. Fig. 6 shows these 
numerical results versus potential width d and Fig. 7 shows those versus Vo/E^y, where E{) is an incident energy and 
Vo is a potential height. It has been shown that, in the opaque case, the Fokker-Planck equation for the distribution 
for the samples can be solved analytically and gives TNeison ~ (= '^wkb) [8], [9]. The parameters adopted in Fig. 6 
give an imaginary wave number k = 1 in the unite of which corresponds to the opaque potential except for very thin 
potential barrier. Thus we can see that TNeison and twkb agree with each other. It is notable that Tvis fits also well 
with them except for thin barrier where the opaqueness condition is broken. The imaginary wave number dependence 
of traversal time is shown in Fig. 7 for a fixed and rather thick potential barrier width. The value of k becomes 
larger than 1 for Vq > 2Eq and in this region TNeison agrees with twkb- On the other hand, in the region Vq < 2Eq, 
Twkb becomes to deviate from TNeison, where the opaqueness condition is not satisfied. However Tyis can reproduce 
the value of TNeison for almost all region. From these two figures, we see, in the opaque case, that TNeison coincide with 
Twkb with respect to its dependence on potential width d and on the imaginary wave number k. While there is an 
obvious reason why the twkb can only applicable to the opaque case, one needs not assume any approximation to 
evaluate TNeison in principle. Therefore the latter may represent an characteristic property of time scale for tunneling 
phenomena not only for the opaque case but also for the translucent case. However, both of these traversal times are 
defined only on the bases of theoretical models, but cannot be checked by experiment so easily. It should be noticed 
that Tvis is connected to the experimental data directly, and the theoretical estimation may be checked by experiment 
rather easily. Thus we think that Tvis can be a good candidate presenting time scale of tunneling phenomena both for 
the opaque case and for the translucent case. 



In this paper, we present a proposal for the estimation of Biittiker-Landauer traversal time based on the visibility 
of transmission current. We analyzed the tunneling phenomena with a time-dependent potential described by Eq. 
(3), and obtained the time-dependent transmission current for a small perturbation V\ and an opaque case. We found 
that the visibility is directly connected to the traversal time, while Biittiker and Landauer proposed that the crossover 
from the low frequency behavior to the high frequency behavior yields the traversal time. Furthermore, this result is 
valid not only for rectangular potential barrier but also for general form of potential to which the WKB approximation 
is applicable. After a brief review of Nelson's quantum mechanics, by which the traversal time is calculated definitely, 
we compared those results with the numerical values obtained from the simulation of Nelson's framework. Both of 
them fit together not only for the opaque case but also for the translucent case and it shows our method is very 
effective to measure experimentally the traversal time. 
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In this Appendix, we recapitulate briefly how to determine the amplitudes of the sidebands at £J ± nfiu) from the 
matching conditions [4], [5], [6]. It is convenient to define the following quantities: 
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APPENDIX A: AMPLITUDES 




(Al) 



(A2) 
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where the prime means a derivative with respect to the coordinate x. At the energy En, we have the matching 
conditions, 



^noe-'"''+A„ e'' 



(A3) 
(A4) 
(A5) 
(A6) 



where a„ = knd/2. Noticing that J„ is proportional to (yi/2?ki;)" and taking only the leading terms, we approximate 
equations (A3), (A4), (A5) and (A6) and obtain the transmission coefficients and reflection coefficients. At the energy 
E, we recover the results of static barrier, the reflection and the transmission coefficients. 



-2(fc^ + K?) sinh nd 



det(A;, k) 



-AikK 
det(fc, k) 



-ikd 



(A7) 
(A8) 



where det(A:, k) is defined by 



det(A;, k) 



{K + ik)e -{n-ik) 
{K-ik)e'^'^ -{n + ik) 

2(/t^ — k^) sinh nd — 4ifc«;cosh Kd. 



(A9) 



Similarly, at the energy £„, we have the transmission and reflection coefficients in the leading order, 

J„(Fi/au;) 2DqS*'-^^)'^/'^ 



and 



Jo(Vi/?iw) det(fc„,K„) 

X {(k^ — fcn/co) sinh K„(i — (k^ — fc„fco)(K„/«;o) sinh/tod 
+ iKn{kn + ko){cosh Kod — coshK„d)} , 

Jo{Vi/hu;) det(fc„,K„) 

X {{k^ — knko) sinh Knd cosh Kod 

— {k^ + knko) (Kn/Ko) cosh Knd sinh Kod 
+iKn{ko — fc„)(l — cosh K„d cosh Korf) 
—i{{koKn/Ko) — A;„/to) sinh K„d sinh kqc^} , 



(AlO) 



(All) 



where det(/u„, Kn) is defined by 

det(A;7i, f^ji) — 



2(/c^ — fc^) sinh Knd — Aiknf^n cosh Knd. 



(A12) 



APPENDIX B: THE VISIBILITY IN A GENERAL POTENTIAL CASE 

We give an expression for the visibility in a general potential case by the use of the WKB approximation. A 
stationary solution ^e{x) with energy E satisfies the Schrodinger equation 
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2m 



+ V{x) 



Stating from the outgoing wave solution in the region III, 



cxp < I 



we have the evanescent wave solution in the region II, 



k{x')dx' 



s 



exp|— «;(a;')da;'| + ^ exp K{x')dx' 



and then, the incoming and reflecting wave solutions 



^\x) = 



exp|-z (^J\{x')dx' -^^"^ 



.4-5^ 



exp|i(^£ ' k{x')dx' - 



where 



S = exp ^— J K{x')dx'^ . 



Using these stationary WKB solutions, we can write down a time dependent solution in the 

For the case of n = 0, we have the solution 



DoJo (^) -j^ 



+ (-^^) exp{i(/r k{x')dx' - f )}] , {x < x^), 

DoJo it) 77hr. f-i «-P{- C 

+ ^ cxp{J^^ K{x')dx'}\ ixi<x< X2), 



/k{x) 

where Sn is the damping factor of the n-th mode, 

Sn = exp 



DoJo it) 77^ expO(/;^ kix')dx' - f )} 



(X2 < X), 



Kn{x)dx 



The coefficients are fixed by the incoming wave normalization. 



-Do Jo 



. 4^0 



For the case of n = 1, we have to consider the two types of wave in the region II, 

(^) Tfe [-i ^-P{- C + f exp{/;^ n{x')dx'} 

+ ^^"^0 (li) vfe [-i e^P{- /x. + f exp{/;^ K,{x')dx'} 



both of which should be matched to the reflecting and the transmitting waves with energy Ei and wave number ki. 
A similar relation holds for n = —1. Requiring the condition that there are no incoming waves in these modes, we 
can determine the coefficients D±i in the region III, 



4 + 52 5±i 



DoJ±i 



TujJ 



4 + 5|i So 

Considering the above results, we get the transmitting wave up to n = ±1, 



4 + 5o'^/fe 



1 f TT 

== exp \ i{ / kQ{x')dx' - - 

•0{X) I Jx2 4 



■iEt/n 



+ 



jQ{Vi/hoj)]] ki{x) [ J^^nkoix') J' 
J_i(Fi/M ' 



ko{x) 

Jo(yi/M V k-i{x) 



exp 



hko{x') 



-dx' \ (1 - E_i)e*' 



where 



4 + ^g S±i 
4 + 5^1 So ■ 



(Bll) 



(B12) 



(B13) 
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FIG. 1. Particles transmitted or reflected at a barrier of height Vb and width d interacting a small modulation Vi cosut can 
absorb or emit modulation quanta hu. The transmitted and reflected waves contain amplitudes at the frequency E/h and the 
sideband frequencies En/Ti. 
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FIG. 2. The transmission probability taking a long time average. 
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FIG. 3. The time dependence of the transmitted currents at a fixed point x = 750(units of l/fco). The potential frequency lj 
is 0.1(units of /cq). Other parameters (static potential height, small modulation amplitude, etc.) are the same values in Fig. 2. 
In this figure, To is the transmission probability in the static potential case, that is, Vi — 0. 
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FIG. 4. Schematical illustration of one-dimensional tunneling in the general potential case. The small modulation Vi cosujt 
exists in the region of II (illustrated with the dashed line) . 
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FIG. 5. Typical transmission sample path calculated by Eq.(33). 





FIG. 6. Comparison of numerical results of traversal times versus potential width d in a rectangular potential barrier. 
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